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Abstract 

A two-dimensional quantum system with anyonic excitations can be considered as a 
quantum computer. Unitary transformations can be performed by moving the excitations 
around each other. Measurements can be performed by joining excitations in pairs and 
observing the result of fusion. Such computation is fault-tolerant by its physical nature. 

A quantum computer can provide fast solution for certain computational problems (e.g. 
factoring and discrete logarithm [0) which require exponential time on an ordinary computer. 
Physical realization of a quantum computer is a big challenge for scientists. One important 
problem is decoherence and systematic errors in unitary transformations which occur in real 
quantum systems. From the purely theoretical point of view, this problem has been solved 
due to Shor's discovery of fault-tolerant quantum computation 0, with subsequent improve- 
ments [§, |4], H |6| . An arbitrary quantum circuit can be simulated using imperfect gates, provided 
these gates are close to the ideal ones up to a constant precision 5. Unfortunately, the threshold 
value of 5 is rather smallQ; it is very difficult to achieve this precision. 

Needless to say, classical computation can be also performed fault-tolerantly. However, it 
is rarely done in practice because classical gates are reliable enough. Why is it possible? Let 
us try to understand the easiest thing — why classical information can be stored reliably on a 
magnetic media. Magnetism arise from spins of individual atoms. Each spin is quite sensitive 
to thermal fluctuations. But the spins interact with each other and tend to be oriented in the 
same direction. If some spin flips to the opposite direction, the interaction forces it to flip back 
to the direction of other spins. This process is quite similar to the standard error correction 
procedure for the repetition code. We may say that errors are being corrected at the physical 
level. Can we propose something similar in the quantum case? Yes, but it is not so simple. 
First of all, we need a quantum code with local stabilizer operators. 

I start with a class of stabilizer quantum codes associated with lattices on the torus and 
other 2-D surfaces || |J. Qubits live on the edges of the lattice whereas the stabilizer operators 
correspond to the vertices and the faces. These operators can be put together to make up a 

1 Actually, the threshold is not known. Estimates vary from 1/300 JtJ to 10~ 6 Q. 
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Hamiltonian with local interaction. (This is a kind of penalty function; violating each stabilizer 
condition costs energy). The ground state of this Hamiltonian coincides with the protected 
space of the code. It is 4 9 -fold degenerate, where g is the genus of the surface. The degeneracy 
is persistent to local perturbation. Under small enough perturbation, the splitting of the ground 
state is estimated as exp(— aL), where L is the smallest dimension of the lattice. This model 
may be considered as a quantum memory, where stability is attained at the physical level rather 
than by an explicit error correction procedure. 

Excitations in this model are anyons, meaning that the global wavefunction acquires some 
phase factor when one excitation moves around the other. One can operate on the ground 
state space by creating an excitation pair, moving one of the excitations around the torus, and 
annihilating it with the other one. Unfortunately, such operations do not form a complete basis. 
It seems this problem can be removed in a more general model (or models) where the Hilbert 
space of a qubit have dimensionality > 2. This model is related to Hopf algebras. 

In the new model, we don't need torus to have degeneracy. An n-particle excited state 
on the plane is already degenerate, unless the particles (excitations) come close to each other. 
These particles are nonabelian anyons, i.e. the degenerate state undergoes a nontrivial unitary 
transformation when one particle moves around the other. Such motion ("braiding") can be 
considered as fault-tolerant quantum computation. A measurement of the final state can be 
performed by joining the particles in pairs and observing the result of fusion. 

Anyons have been studied extensively in the field-theoretic context || [H], [11], [12], |I3| . So, 
I hardly discover any new about their algebraic properties. However, my approach differs in 
several respects: 

• The model Hamiltonians are different. 

• We allow a generic (but weak enough) perturbation which removes any symmetry of the 
Hamiltonian.^ 

• The language of ribbon and local operators (see Sec. |5.2j ) provides unified description of 
anyonic excitations and long range entanglement in the ground state. 

An attempt to use one- dimensional anyons for quantum computation was made by G. Castag- 



noli and M. Rasetti [14], but the question of fault-tolerance was not considered. 



1 Toric codes and the corresponding Hamiltonians 

Consider a k x k square lattice on the torus (see fig. [I]). Let us attach a spin, or qubit, to 
each edge of the lattice. (Thus, there are n = 2k 2 qubits). For each vertex s and each face p, 
consider operators of the following form 

a s = n °* b v= n 

j'6star(s) j€boundary(p) 

These operators commute with each other because star(s) and boundary(p) have either or 2 
common edges. The operators A s and B p are Hermitian and have eigenvalues 1 and —1. 

2 Some local symmetry still can be established by adding unphysical degrees of freedom, see Sec. ||. 
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Figure 1: Square lattice on the torus 



Let M be the Hilbert space of all n = 2k 2 qubits. Define a protected subspace £ C N as 
follows^ 

£ = { |0 G AT : A s \0 = |0, SpIO = 10 for all S ,p } (2) 

This construction gives us a definition of a quantum code TOR(fc), called a tone code |], ||. 
The operators A s , B p are the stabilizer operators of this code. 

To find the dimensionality of the subspace £, we can observe that there are two relations 
between the stabilizer operators, EL^s = 1 an d Yl p B P = 1- So, there are m = 2k 2 — 
2 independent stabilizer operators. It follows from the general theory of additive quantum 
codes |Lj|, [T6| that dim£ = 2 n ~ m = 4. 

However, there is a more instructive way of computing dim£. Let us find the algebra L(£) 
of all linear operators on the space £ — this will give us full information about this space. 
Let T C L(yV) be the algebra of operators generated by A s , B p . Clearly, L(£) = Q/I, where 
Q ~D T is the algebra of all operators which commute with A s , B p , and I C Q is the ideal 
generated by A s — 1, B p — 1. The algebra £ is generated by operators of the form 

z = n°i x = n°i 

where c is a loop (closed path) on the lattice, whereas d is a cut, i.e. a loop on the dual lattice 
(see fig. H). If a loop (or a cut) is contractible then the operator Z is a product of 5 p , hence 
Z = 1 (mod X). Thus, only non-contractible loops or cuts are interesting. It follows that the 
algebra L(£) is generated by 4 operators Z ly Z 2 , Xx, X 2 corresponding to the loops c z i, c z2 
and the cuts c x i, c x2 (see fig. H). The operators Z ly Z 2 , Xx, X 2 have the same commutation 
relations as a\, a 2 , af, a 2 . We see that each quantum state |0 G £ corresponds to a state of 
2 qubits. Hence, the protected subspace £ is 4- dimensional. 

In a more abstract language, the algebra T corresponds to 2-boundaries and O-coboundaries 
(with coefficients from r L 2 ),Q corresponds to 1-cycles and 1-cocycles, and L(£) corresponds to 
1-homologies and 1-cohomologies. 

There is also an explicit description of the protected subspace which may be not so useful 
but is easier to grasp. Let us choose basis vectors in the Hilbert space M by assigning a label 
Zj = 0, 1 to each edge j. Q The constraints B p \£) = |0 say that the sum of the labels at the 
boundary of a face should be zero (mod 2). More exactly, only such basis vectors contribute to 

3 We will show that this subspace is really protected from certain errors. Vectors of this subspaces are 
supposed to represent "quantum information" , like codewords of a classical code represent classical information. 
means "spin up", 1 means "spin down". The Pauli operators cr 2 , a x have the standard form in this basis. 
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Figure 2: Loops on the lattice and the dual lattice 



a vector from the protected subspace. Such a basis vector is characterized by two topological 
numbers: the sums of Zj along the loops c z \ and c z2 . The constraints A p \£) = |£) say that all 
basis vectors with the same topological numbers enter |£) with equal coefficients. Thus, for 
each of the 4 possible combinations of the topological numbers v i, t>2, there is one vector from 
the protected subspace, 

Of course, one can also create linear combinations of these vectors. 

Now we are to show that the code TOR(fc) detects k — 1 errors^] (hence, it corrects [_^rj 
errors). Consider a multiple error 

E = a( ai ,..., a*; = IlK) Qj II («i,/9i = 0,l) 

This error can not be detected by syndrome measurement (i.e. by measuring the eigenvalues of 
all A 8 , Bj) if and only if E G Q. However, if E G T then E\£) = |^) for every |^) G £. Such an 
error is not an error at all — it does not affect the protected subspace. The bad case is when 
E G Q but E G" J- ' . Hence, the support of E should contain a non-contractible loop or cut. It 
is only possible if | Supp(£')| > k. (Here Supp(i?) is the set of j for which ay ^ or (3j ^ 0). 

One may say that the toric codes have quite poor parameters. Well, they are not "good" 
codes in the sense of [17]. However, the code TOR(fc) corrects almost any multiple error of size 
0{k 2 ). (The constant factor in 0(. . .) is related to the percolation problem). So, these codes 
work if the error rate is constant but smaller than some threshold value. The nicest property 
of the codes TOR(fc) is that they are local check codes. Namely, 

1. Each stabilizer operators involves bounded number of qubits (at most 4). 

2. Each qubit is involved in a bounded number of stabilizer operators (at most 4). 



In the theory of quantum codes, the word "error" is used in a somewhat confusing manner. Here it means 
a single qubit error. In most other cases, like in the formula below, it means a multiple error, i.e. an arbitrary 
operator E e L(J\f). 
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3. There is no limit for the number of errors that can be corrected. 

Also, at a constant error rate, the unrecoverable error probability goes to zero as exp(— ak). 

It has been already mentioned that error detection involves syndrome measurement. To 
correct the error, one needs to find its characteristic vector (a>i, . . . , a n ] Pi, ... , (3 n ) out of the 
syndrome. This is the usual error correction scheme. A new suggestion is to perform error 
correction at the physical level. Consider the Hamiltonian 

h = - J2^ s - J2 b p ( 4 ) 

s p 

Diagonalizing this Hamiltonian is an easy job because the operators A s , B p commute. In 
particular, the ground state coincides with the protected subspace of the code TOR(&); it is 
4- fold degenerate. All excited states are separated by an energy gap AE > 2, because the 
difference between the eigenvalues of A s or B p equals 2. This Hamiltonian is more or less 
realistic because in involves only local interactions. We can expect that "errors", i.e. noise- 
induced excitations will be removed automatically by some relaxation processes. Of course, 
this requires cooling, i.e. some coupling to a thermal bath with low temperature (in addition 
to the Hamiltonian (|])). 

Now let us see whether this model is stable to perturbation. (If not, there is no practical 
use of it). For example, consider a perturbation of the form 

j j<p 

It is important that the perturbation is local, i.e. each term of it contains a small number of 
a (at most 2). Let us estimate the energy splitting between two orthogonal ground states of 
the original Hamiltonian, |£) e £ and G C. We can use the usual perturbation theory 
because the energy spectrum has a gap. In the m-th order of the perturbation theory, the 
splitting is proportional to (^\V m \r]) or (^IV 71 )^) — (r]\V m \ri). However, both quantities are 
zero unless V m contains a product of cr| or crj along a non-contractible loop or cut. Hence, the 
splitting appears only in the |~fc/2~|-th or higher orders. As far as all things (like the number 
of the relevant terms in V m ) scale correctly to the thermodynamic limit, the splitting vanishes 
as exp(— ak). A simple physical interpretation of this result is given in the next section. (Of 
course, the perturbation should be small enough, or else a phase transition may occur). 

Note that our construction is not restricted to square lattices. We can consider an arbitrary 
irregular lattice, like in fig. ^|. Moreover, such a lattice can be drawn on an arbitrary 2-D 
surface. On a compact orientable surface of genus g, the ground state is 4 9 -fold degenerate. In 
this case, the splitting of the ground state is estimated as exp(— aL), where L is the smallest 
dimension of the lattice. We see that the ground state degeneracy depends on the surface 
topology, so we deal with topological quantum order. On the other hand, there is a finite 
energy gap between the ground state and excited states, so all spatial correlation functions 
decay exponentially. This looks like a paradox — how do parts of a macroscopic system know 
about the topology if all correlations are already lost at small scales? The answer is that there 
is long-range entanglement | which can not be expressed by simple correlation functions like 
(cr^of). This entanglement reveals itself in the excitation properties we are going to discuss. 

6 Entanglement is a special, purely quantum form of correlation. 
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Figure 3: Strings and particles 

2 Abelian anyons 

Let us classify low-energy excitations of the Hamiltonian (^). An eigenvector of this Hamiltonian 
is an eigenvector of all the operators A s , B p . An elementary excitation, or particle occurs if only 
one of the constraints A s \£) = B p £) = |£) is violated. Because of the relations Y[ s A s = 1 
and rip-^p = 1) ^ * s impossible to create a single particle. However, it is possible to create a 
two-particle state of the form \ip z {t)) = S z (t)\£) or \ip x (t')) = S x (t')\£), where |£) is an arbitrary 
ground state, and 

si® = = ( 5 ) 

jet jet' 

(see fig. |3|). In the first case, two particles are created at the endpoints of the "string" (non- 
closed path) t. Such particles live on the vertices of the lattice. We will call them z-type 
particles, or "electric charges". Correspondingly, x-type particles, or "magnetic vortices" live 
on the faces. The operators S z (t), S x (t') are called string operators. Their characteristic 
property is as follows: they commute with every A s and B p , except for few ones (namely, 2) 
corresponding to the endpoints of the string. Note that the state \ip z (t)) = S z (t)\£) depends 
only on the homotopy class of the path t while the operator S z (t) depends on t itself. 

Any configuration of an even number of z-type particles and an even number of x-type 
particles is allowed. We can connect them by strings in an arbitrary way. Each particle 
configuration defines a 4 5 -dimensional subspace in the global Hilbert space M. This subspace 
is independent of the strings but a particular vector S ai (ti) ■ ■ ■ S am (t m ) |£) depends on ti, . . . , t m . 
If we draw these strings in a topologically different way, we get another vector in the same 4 5 - 
dimensional subspace. Thus, the strings are unphysical but we can not get rid of them in our 
formalism. 

Let us see what happens if these particles move around the torus (or other surface). Moving 
a z-type particle along the path c z \ or c Z 2 (see fig. ||) is equivalent to applying the operator Z\ 
or Z 2 - Thus, we can operate on the ground state space by creating a particle pair, moving one 
of the particles around the torus, and annihilating it with the other one. Thus we can realize 
some quantum gates. Unfortunately, too simple ones — we can only apply the operators a z 
and a x to each of the 2 (or 2g) qubits encoded in the ground state. 

Now we can give the promised physical interpretation of the ground state splitting. In the 
presence of perturbation, the two-particle state \ip z (t)) is not an eigenstate any more. More 
exactly, both particles will propagate rather than stay at the same positions. The propagation 
process is described by the Schrodinger equation with some effective mass m z . (x-type particles 
have another mass m x ). In the non-perturbed model, m z = m x = oo. There are no real particles 
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Figure 4: An x-type particle moving around a z-type particle 



in the ground state, but they can be created and annihilate virtually. A virtual particle can 
tunnel around the torus before annihilating with the other one. Such processes contribute terms 
b z \Zi, b z2 Z 2 , b x \Xi, b x2 X 2 to the ground state effective Hamiltonian. Here b a i ~ exp(— a a Li) is 
the tunneling amplitude whereas a a ~ yhnKE is the imaginary wave vector of the tunneling 
particle. 

Next question: what happens if we move particles around each other? (For this, we don't 
need a torus; we can work on the plane). For example, let us move an x-type particle around 
a z-type particle (see fig. f|). Then 

^initial) = S Z {t)\r{q)) , = S X (c)S Z (t)\r(q)) = "^initial) 

because S x (c) and S z (t) anti-commute, and S x (c)\ip x (q)) = \ip x (q)). We see that the global 
wave function (= the state of the entire system) acquires the phase factor —1. It is quite 
unlike usual particles, bosons and fermions, which do not change their phase in such a process. 
Particles with this unusual property are called abelian anyons. More generally, abelian anyons 
are particles which realize nontrivial one-dimensional representations of (colored) braid groups. 
In our case, the phase change can be also interpreted as an Aharonov-Bohm effect. It does not 
occur if both particles are of the same type. 

Note that abelian anyons exist in real solid state systems, namely, they are intrinsicly re- 
lated to the fractional quantum Hall effect [13]. However, these anyons have different braiding 
properties. In the fractional quantum Hall system with filling factor p/q, there is only one 
basic type of anyonic particles with (real) electric charge 1/q. (Other particles are thought to 
be composed from these ones). When one particle moves around the other, the wave function 
acquires a phase factor exp(2iri/q). 

Clearly, the existence of anyons and the ground state degeneracy have the same nature. They 
both are manifestations of a topological quantum order, a hidden long-range order that can not 
be described by any local order parameter. (The existence of a local order parameter contradicts 
the nature of a quantum code — if the ground state is accessible to local measurements then 
it is not protected from local errors). It seems that the anyons are more fundamental and can 
be used as a universal probe for this hidden order. Indeed, the ground state degeneracy on the 
torus follows from the existence of anyons ||19|| . Here is the original Einarsson's proof applied 



to our two types of particles. 

We derived the ground state degeneracy from the commutation relations between the op- 
erators Zi, Z 2 , Xi, X 2 . These operators can be realized by moving particles along the loops 
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Figure 5: A fly-over crossing geometry for a 2-D electron layer 



Czii Cz2i Cxi) c x2- These loops only exist on the torus, not on the plane. Consider, however, 
the process in which an x-type particle and a £-type particle go around the torus and then 
trace their paths backward. This corresponds to the operator W = X^ 1 Z^~ 1 XiZ 1 which can be 
realized on the plane. Indeed, we can deform particles' trajectories so that one particle stays 
at rest while the other going around it. Due to the anyonic nature of the particles, W = — 1. 
We see that X\ and Z\ anti-commute. 

The above argument is also applicable to the fractional quantum Hall anyons fl9|| . The 
ground state on the torus is g-fold degenerate, up to the precision ~ exp(— L/l ), where l 
is the magnetic length. This result does not rely on the magnetic translational symmetry or 
any other symmetry. Rather, it relies on the existence of the energy gap in the spectrum 
(otherwise the degeneracy would be unstable to perturbation). Note that holes (=punctures) 
in the torus do not remove the degeneracy unless they break the nontrivial loops c x i, c X 2, c z ±, 
c Z 2- The fly-over crossing geometry (see fig. |) is topologically equivalent to a torus with 2 
holes, but it is almost flat. In principle, such structure can be manufactured f\, cooled down 
and placed into a perpendicular magnetic field. This will be a sort of quantum memory — it 
will store a quantum state forever, provided all anyonic excitation are frozen out or localized. 
Unfortunately, I do not know any way this quantum information can get in or out. Too few 
things can be done by moving abelian anyons. All other imaginable ways of accessing the 
ground state are uncontrollable. 

3 Materialized symmetry: is that a miracle? 

Anyons have been studied extensively in the gauge field theory context || [10, 0, |T^] . However, 
we start with quite different assumptions about the Hamiltonian. A gauge theory implies a 
gauge symmetry which can not be removed by external perturbation. To the contrary, our model 
is stable to arbitrary local perturbations. It is useful to give a field-theoretic interpretation 
of this model. The edge labels Zj (measurable by crj) correspond to a Z 2 vector potential, 
whereas crj corresponds to the electric field. The operators A s are local gauge transformations 
whereas B p is the magnetic field on the face p. The constraints A s \£) = |£) mean that the 
state |£) is gauge- invariant. Violating the gauge invariance is energetically unfavorable but not 
forbidden. The Hamiltonian (which includes Hq and some perturbation V) need not obey the 
gauge symmetry. The constraints B p \^) = |£) mean that the gauge field corresponds to a flat 
connection. These constraints are not strict either. 

7 It is not easy. How will the two layers (the two crossing "roads", one above the other) join in a single 
crystal layout? 
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Despite the absence of symmetry in the Hamiltonian H = H + V, our system exhibits two 
conservation laws: electric charge and magnetic charge (i.e. the number of vortices) are both 
conserved modulo 2. In the usual electrodynamics, conservation of electric charge is related to 
the local (=gauge) U(l) symmetry. In our case, it should be a local Z 2 symmetry for electric 
charges and another Z 2 symmetry for magnetic vortices. So, our system exhibits a dynamically 
created Z 2 x Z 2 symmetry which appears only at large distances where individual excitations 
are well-defined. 

Probably, the reader is not satisfied with this interpretation. Really, it creates a new puzzle 
rather than solve an old one. What is this mysterious symmetry? How do symmetry operators 
look like at the microscopic level? The answer sounds as nonsense but it is true. This symmetry 
(as well as any other local symmetry) can be found in any Hamiltonian if we introduce some 
unphysical degrees of freedom. So, the symmetry is not actually being created. Rather, an 
artificial symmetry becomes a real one. 

The new degrees of freedom are spin variables v s ,w p — 0, 1 for each vertex s and each face 

p. The vertex spins will stay in the state 2 -1 / 2 ^|0) + whereas the face spins will stay in 
the state |0). So, all the extra spins together stay in a unique quantum state \(). Obviously, 
a s 10 = 10 an d CpIO = 10) f° r ever Y vertex s and every face p. From the mathematical point 
of view, we have simply defined an embedding of the space N into a larger Hilbert space T of 
all the spins, \ip) i— > \ip) ® \Q. So we may write N C T. We will call N the physical space (or 
subspace), T the extended space. Physical states (i.e. vectors G N) are characterized by 
the equations 

cm = m , o* P w) = iv) 

for every vertex s and face p. 

Now let us apply a certain unitary transformation U on the extended space T. This trans- 
formation is just a change of the spin variables, namely 

v s h-> v s , Zj h-> Zj + ^2 v s , w p l-> w p + ^2 z i 

s=endpoint( i j) boundary (p) 

(all sums are taken modulo 2). The physical subspace becomes N' = UN. Vectors \ip) G N 
are invariant under the following symmetry operators 

P s = Ua x s tf = a*A s Q p = Ua z p U^ = a z p B p (7) 

The transformed Hamiltonian H' = U HU^ commutes with these operators. It is defined up to 
the equivalence P s = 1, Q p = 1. In particular, 

H' = UH rf = H' Q = - J2°s -J2 a P ( 8 ) 

s p 

In the field theory language, the vertex variables v s (or the operators af) are a Higgs field. 
The operators P s are local gauge transformations. Thus, an arbitrary Hamiltonian can be 
written in a gauge-invariant form if we introduce additional Higgs fields. Of course, it is a 
very simple observation. The real problem is to understand how the artificial gauge symmetry 
"materialize", i.e. give rise to a physical conservation law. 
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Electric charge at a vertex s is given by the operator a x s . The total electric charge on a 
compact surface is zeroQ because fls^s = This is not a physically meaningful statement 
as it is. It is only meaningful if there are discrete charged particles. Then the charge is also 
conserved locally, in every scattering or fusion process. It is difficult to formulate this property 
in a mathematical language, but, hopefully, it is possible. (The problem is that particles are 
generally smeared and can propagate. Physically, particles are well-defined if they are stable and 
have finite energy gap). Alternatively, one can use various local and nonlocal order parameters 
to distinguish between phases with an unbroken symmetry, broken symmetry or confinement. 

The artificial gauge symmetry materialize for the Hamiltonian (J8|) but this is not the case 
for every Hamiltonian. Let us try to describe possible symmetry breaking mechanisms in terms 
of local order parameters. If the gauge symmetry is broken then there is a nonvanishing vacuum 
average of the Higgs field, <fi(s) = [af) 7^ 0. Electric charge is not conserved any more. In other 
words, there is a Bose condensate of charged particles. Although the second Z 2 symmetry is 
formally unbroken, free magnetic vortices do not exist. More exactly, magnetic vortices are 
confined. (The duality between symmetry breaking and confinement is well known P5|). It is 
also possible that the second symmetry is broken, then electric charges are confined. From the 
physical point of view, these two possibilities are equivalent: there is no conservation law in 
the system.^ 

An interesting question is whether magnetic vortices can be confined without the gauge 
symmetry being broken. Apparently, the answer is "no" . The consequence is significant: electric 
charges and magnetic vortices can not exist without each other. It seems that materialized 
symmetry needs better understanding; as presented here, it looks more like a miracle. 



4 The model based on a group algebra 

From now on, we are constructing and studying nonabelian anyons which will allow universal 
quantum computation. 

Let G be a finite (generally, nonabelian) group. Denote by Ji = C[G] the corresponding 
group algebra, i.e. the space of formal linear combinations of group elements with complex 

coefficients. We can consider 7i as a Hilbert space with a standard orthonormal basis : 

g G Gj. The dimensionality of this space is iV = We will work with "spins" (or "qubits") 

taking values in this space.0 Remark: This model can be generalized. One can take for 7i 
any finite-dimensional Hopf algebra equipped with a Hermitian scalar product with certain 
properties. However, I do not want to make things too complicated. 

To describe the model, we need to define 4 types of linear operators, L 9 + , L 9 _, T+, acting 
on the space H. Within each type, they are indexed by group elements, g G G or h G G. They 
act as follows 

L 9 + \z) = \gz) T$\z) = 8 hjZ \z) 

Ll\z) = \zg- 1 ) Tl\z) = 8 h -r z \z) (9) 

8 Strictly speaking, the electric charge is not a numeric quantity; rather, it is an irreducible representation 
of the group Z2. "Zero" refers to the identity representation. 

9 The two possibilities only differ if an already materialized symmetry breaks down at much large distances 
(lower energies). 

10 In the field theory language, the value of a spin can be interpreted as a G gauge field. However, we do not 
perform symmetrization over gauge transformations. 
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Figure 6: Generic lattice and the orientation rules for the operators L 9 ± and 

(In the Hopf algebra context, the operators L g + , L g _, T+, correspond to the left and right 
multiplications and left and right comultiplication, respectively). These operators satisfy the 
following commutation relations 

I 9 T h = T gh L g L g T h = T hg l L g 

+ + + + + - - + (1Q) 

If T h = T hg L 9 L g T h = T gh L g 

Now consider an arbitrary lattice on an arbitrary orientable 2-D surface, see fig. |6|. (We 
will mostly work with a plane or a sphere, not higher genus surfaces). Corresponding to each 
edge is a spin which takes values in the space Ti. Arrows in fig. || mean that we choose some 
orientation for each edge of the lattice. (Changing the direction of a particular arrow will be 
equivalent to the basis change \z) i— »■ \z^) for the corresponding qubit). Let j be an edge of the 
lattice, s one of its endpoints. Define an operator L 9 (j, s) = L±(j) as follows. If s is the origin 
of the arrow j then L 9 (j,s) is L g _(j) (i.e. L g _ acting on the j-th spin), otherwise it is L g + (j). 
This rule is represented by the diagram at the right side of fig. |6|. Similarly, if p is the left (the 
right) ajacent face of the edge j then T h (j,p) is (resp. T+) acting on the j-th spin. 

Using these notations, we can define local gauge transformations and magnetic charge op- 
erators corresponding to a vertex s and an adjacent face p (see fig. ||). Put 

A g (s,p) = A g (s) = n L 9 (j,s) 

j'Gstar(s) 

(11) 

B h (s,p) = E 11 T h -( Jm ,p) 

hi---h^=h m=l 

where ji, ■ ■ ■ ,jk are the boundary edges of p listed in the counterclockwise order, starting from, 
and ending at, the vertex s. (The sum is taken over all combinations of hi, . . . , hk G G, such that 
hi ■ ■ ■ hfc = h. Order is important here!). Although A g (s,p) does not depend onp, we retain this 
parameter to emphasize the duality between A g (s,p) and B h (s,p). These operators generate 
an algebra T> = D(G), Drinfield's quantum double [2D| of the group algebra C[G]. It will play 



a very important role below. Now we only need two symmetric combinations of A g (s,p) and 
Bh(s,p), namely 

A(s) = N-'Y^M^P) B{p) = Bi(s,p) (12) 

geG 
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In the Hopf algebra setting, A g (s,p) does depend on p. 
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where N = \G\. Both A(s) and B(p) are projection operators. (A(s) projects out the states 
which are gauge invariant at s, whereas B(p) projects out the states with vanishing magnetic 
charge at p). The operators A(s) and B(p) commute with each other.0 Also A(s) commutes 
with A(s'), and B(p) commutes with B(p') for different vertices and faces. In the case G = Z 2 , 
these operators are almost the same as the operators ([!]), namely A(s) = \{A S + 1), B{p) = 

|(s P + i).Q 

At this point, we have only defined the global Hilbert space M (the tensor product of many 
copies of 7i) and some operators on it. Now let us define the Hamiltonian. 

Ho = YfX-Ms)) + £(1-B(p)) (13) 

s p 

It is quite similar to the Hamiltonian (Q). As in that case, the space of ground states is given 
by the formula 

C = { |0 G N : A(s)\0 = |0, B(p)\0 = |0 for all s,p} (14) 

The corresponding energy is 0; all excited states have energies > 1. 

It is easy to work out an explicit representation of ground states similar to eq. ([|). The 
ground states correspond 1-to-l to flat G- connect ions, defined up to conjugation, or super- 
positions of those. So, the ground state on a sphere is not degenerate. However, particles 
(excitations) have quite interesting properties even on the sphere or on the plane. (We treat 
the plane as an infinitely large sphere). The reader probably wants to know the answer first, 
and then follow formal calculations. So, I give a brief abstract description of these particles. It 
is a mixture of general arguments and details which require verification. 

The particles live on vertices or faces, or both; in general, one particle occupies a vertex 
and an adjacent face same time. A combination of a vertex and an adjacent face will be called 
a site. Sites are represented by dotted lines in fig. [7]. (The dashed lines are edges of the dual 
lattice). 

Consider n particles on the sphere pinned to particular sites Xi, . . . ,x n at large distances 
from each other. The space C[n] = C(xi, . . . , x n ) of n-particle states has dimensionality JV 2 ^ -1 ', 
including the ground stateQ Not all these states have the same energy. Even more splitting 
occurs under perturbation, but some degeneracy still survive. Of course, we assume that the 
perturbation is local, i.e. it can be represented by a sum of operators each of which acts only on 
few spins. To find the residual degeneracy, we will study the action of such local operators on 
the space C[n). Local operators generate a subalgebra V[n] C L(£[n]). Elements of its center, 
C[n], are conserved classical quantities; they can be measured once and never change. (More 
exactly, they can not be changed by local operators). As these classical variables are locally 
measurable, we interpret them as particle's types. It turns out that the types correspond 1-to-l 
to irreducible representations of the algebra V, the quantum double. Thus, each particle can 
belong to one of these types. The space C[n] and the algebra V[n] split accordingly: 

C[n] = © V[n] = t (v dl ,..., dn QUU,..., dn )) (15) 

di,...,d n d\,...,d n 

12 This is not obvious. Use the commutation relations (|l(]) to verify this statement. 

13 Here A s and B p are the notations from Sec. [IJ we will not use them any more. 

14 The absence of particle at a given site is regarded as a particle of special type. 
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where d m is the type of the m-th particle. The "classical" subalgebra C[n] is generated by the 
projectors onto C dl) ... jdn . 

But this is not the whole story The subspace C dl ,... t d n splits under local perturbations from 
'Pd 1 ,...,d„- By a general mathematical argument,^ this algebra can be characterized as follows 

C dl ,...,d n = ^dl,...,dn ® M dl> ..., dn 'Pd 1 ,...,dn = H^d U ...,dn) ( 16 ) 

The space KL dlr __ >dn corresponds to local degrees of freedom. They can be defined independently 
for each particle. So, K. dlj ___ )dn = JC dl (S> • - • ® /Q n , where JC dm is the space of "subtypes" 
(internal states) of the m-th particle. Like the type, the subtype of a particle is accessible by 
local measurements. However, it can be changed, while the type can not. 

The most interesting thing is the protected space M.d 1 ,...,d n - It is not accessible by local 
measurements and is not sensitive to local perturbations, unless the particles come close to each 
other. This is an ideal place to store quantum information and operate with it. Unfortunately, 
the protected space does not have tensor product structure. However, it can be described as 
follows. Associated with each particle type a is an irreducible representation U d of the quantum 
double V. Consider the product representation U dl <8> ■ ■ • <8> U dn and split it into components 
corresponding to different irreducible representations. The protected space is the component 
corresponding to the identity representation. 

If we swap two particles or move one around the other, the protected space undergoes some 
unitary transformation. Thus, the particles realize some multi-dimensional representation of the 
braid group. Such particles are called nonabelian anyons. Note that braiding does not affect the 
local degrees of freedom. If two particles fuse, they can annihilate or become another particle. 
The protected space becomes smaller but some classical information comes out, namely, the 
type of the new particle. So, the we can do measurements on the protected space. Finally, if we 
create a new pair of particles of definite types, it always comes in a particular quantum state. 
So, we have a standard toolkit for quantum computation (new states, unitary transformations 
and measurements), except that the Hilbert space does not have tensor product structure. 
Universality of this toolkit is a separate problem, see Sec. |7|. 

Our model gives rise to the same braiding and fusion rules as gauge field theory models [TO, 



nfl . The existence of local degrees of freedom (subtypes) is a new feature. These degrees of 



freedom appear because there is no explicit gauge symmetry in our model. 



5 Algebraic structure 

5.1 Particles and local operators 

This subsection is also rather abstract but the claims we do are concrete. They will be proven 



in Sec. 



As mentioned above, the ground state of the Hamiltonian (|13l) is not degenerate (on the 
sphere or on the plane regarded as an infinitely large sphere). Excited states are characterized 
by their energies. The energy of an eigenstate \ip) is equal to the number of constraints (A(s) — 
l)\ip) = or (B(p) — = which are violated. Complete classification of excited states is 
a difficult problem. Instead of that, we will try to classify elementary excitations, or particles. 

15 Vd u ...,d n is a subalgebra of L(£<2i,...,d„) with a trivial center, closed under Hermitian conjugation. 
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Let us formulate the problem more precisely. Consider a few excited "spots" separated by 
large distances. Each spot is a small region where some of the constrains are violated. The 
energy of a spot can be decreased by local operators but, generally, the spot can not disappear. 
Rather, it shrinks to some minimal excitation (which need not be unique). We will see (in 
Sec. |5.4|) that any excited spot can be transformed into an excitation which violates at most 
2 constraints, A(s) — 1 = and B(p) — 1 = 0, where s is an arbitrary vertex, and p is an 
adjacent face. Such excitations are be called elementary excitations, or particles. Note that 
definition of elementary excitations is a matter of choice. We could decide that an elementary 
excitation violates 3 constraints. Even with our definition, the "space of elementary excitations" 
is redundant. 

By the way, the space of elementary excitations is not well defined because such an excitation 
does not exist alone. More exactly, the only one-particle state on the sphere is the ground state. 
(This can be proven easily). The right thing is the space of two-particle excitations, C(a,b). 
Here a = (s,p) and b = (s',p') are the sites occupied by the particles. (Recall that a site is 
a combination of a vertex and an adjacent face). The projector onto C(a,b) can be written 
as Ylryts s' ^( r ) Tli^pp' B(l)- Note that introducing a third particle (say, c) will not give more 
freedom for any of the two. Indeed, b and c can fuse without any effect on a. 

Let us see how local operators act on the space C(a, b). In this context, a local operator is 
an operator which acts only on spins near a (or near b). Besides that, it should preserve the 
subspace C(a, b) C J\f and its orthogonal complement. (A/* is the space of all quantum states). 
Example: the operators A g (a) and Bh(a), where a = (s,p), commute with A(r), B(l) for all 
r ^ s and / ^ p. Hence, they commute with the projector onto the subspace C(a,b). These 
operators generate an algebra T>(a) C L(AA). It will be shown in Sec. [5.4| that T>(a) includes 
all local operators acting on the space C(a, b), and the action of T>(a) on C(a,b) is exact (i.e. 
different operators act differently). 

Actually, the algebra T>(a) = V does not depend on a, only the embedding V — ► L(A/") does. 
This algebra is called the quantum double of the group G and denoted by D(G). Its structure 
is determined by the following relations between the operators A g = A g (a) and Bh = Bh{a) 

Af Ag = Afg B h Bi = Sh,iB h AgB h = Bg h g-lAg (17) 



The operators D^ >g ) = BhA g form a linear basis of V. (In [|K], [D] these operators were denoted 

by h \ ). The following multiplication rules hold 

g 

D (h U9l )D( h2 ,g 2 ) = S hugih2g -1 D {hugig2) 

This identity can be also written in a symbolic tensor form, with h and g being combined into 
one index: 

D m D n = £lm n D k ^(htgi) (fta^,) = S h 1 ,g 1 h 2 9i 1 S 9,9i92 ( 18 ) 

(summation over k is implied). Actually, T> is not only an algebra, it is a quasi-triangular Hopf 
algebra, see Sees. |5TB| . 

Note that V = V{a) is closed under Hermitian conjugation (in L(A/")) which acts as follows 

4 = V 1 B{ = B h D\ Kg) = D ig -^ g>g ^ (19) 
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Thus, T> = T>(a) is a finite-dimensional C*-algebra. Hence it has the following structure 



V = 0L(/Q) (20) 

d 

where d runs over all irreducible representations of T>. We can interpret d as particle's type.0 
The absence of particle corresponds to a certain one-dimensional representation called the 
identity representation. More exactly, the operators D^, g ) act on the ground state |£) as follows 

Dk\0 = £k\0 where £ {Kg) = S hjl (21) 

The "space of subtypes" , /Q actually characterize the redundancy of our definition of el- 
ementary excitations. However, this redundancy is necessary to have a nice theory of ribbon 



operators (see Sec. |5.2|) . 



Irreducible representations of T> can be described as follows ||10|| . Let u G G be an arbitrary 
element, C = {gug^ 1 : g G G} its conjugacy class, E = {g G G : gu = ug} its centralizer. 
There is one irreducible representation d = (C, x) for each conjugacy class C and each irre- 
ducible representation x of the group E (see below). It does not matter which element u G C 
is used to define E. The conjugacy class C can be interpreted as magnetic charge whereas x 
corresponds to electric charge. For example, consider the group S3 (the permutation group of 
order 3). It has 3 conjugacy classes of order 1, 2 and 3, respectively. So, the algebra D(S , 3) has 
irreducible representations of dimensionalities 1,1,2; 2,2,2; 3,3. 

The simplest case is when x is the identity representation, i.e. the particle carries only 
magnetic charge but no electric charge. Then the subtypes can be identified with the elements 
of C, i.e. the corresponding space (denoted by Be) has a basis {\v) : v G C}. The local 
operators act on this space as follows 

D(h,g)\v) = 6 h ,gvg-i \gvg~ 1 ) (22) 

Now consider the general case. Denote by Wf = the irreducible action of / G E on an 
appropriate space A x . Choose arbitrary elements q v G G such that q v uq~ l = v for each v G C. 
Then any element g G G can be uniquely represented in the form g = q v f, where v G C and 
/ G E. We can define a unique action of V on Be <8> A x , such that 

Bh(\v)®\7i)) = 8 htV \v)®\ri) (23) 
4*/(l«>®to>) = \v)®W f \r 1 ) (veC,feE) 

More generally, D [h)g) (\v) g> \tj)J = 5 h;gvg -i \gvg~ 1 ) <g> W f \r]), where / = q v {q gvg ^)~ l g. This 
action is irreducible. 



5.2 Ribbon operators 

The next task is to construct a set of operators which can create an arbitrary two-particle state 
from the ground state. I do not know how to deduce an expression for such operators; I will 

16 Caution. The local operators should not be interpreted as symmetry transformations. The true symmetry 
transformations, so-called topological operators, will be defined in Sec. |(| Mathematically, they are described 
by the same algebra T>, but their action on physical states is different. 
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Figure 7: A ribbon on the lattice 



just give an answer and explain why it is correct. In the abelian case (see Sec. |2|) there were 
two types of such operators which corresponded to paths on the lattice and the dual lattice, 
respectively. In the nonabelian case, we have to consider both types of paths together. Thus, 
the operators creating a particle pair are associated with a ribbon (see fig. ^). The ribbon 
connects two sites at which the particles will appear (say, a = (s,p) and b = (s',p')). The 
corresponding operators act on the edges which constitute one side of the ribbon (solid line), 
as well as the edges intersected by the other side (dashed line). 

For a given ribbon t, there are N 2 ribbon operators F^ h ' 9 \t) indexed by g, h G G. They act 
as follows^ 



X'2 



X 3 



yi 



V3 



(24) 



X\ 



X-2 



X 3 



J g, X1X2X3 



hyi 



x x 1 hxi y 2 



>ix 2 ) 1 h(xix 2 )y 3 



These operators commute with every projector A(r), B(l), except for r = s,s' and I = p,p'. 
This is the first important property of ribbon operators. 

The operators F^ h ' 9 \t) depend on the ribbon t. However, their action on the space C(a,b) 

17 Horizontal and vertical arrows are the two types of edges. Each of the two diagrams (6 arrows with labels) 
stand for a particular basis vector 
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depends only on the topological class of the ribbon This is also true for a multi-particle excitation 
space C(x\, . . . , x n ). More exactly, consider two ribbons, t and q, connecting the sites X\ — a 
and x 2 = b The actions of F^ h ' 9 \t) and F^ h,9 \q) on C(xi, . . . , x n ) coincide provided none of the 




sites X3, . . . ,x n lie on or between the ribbons. This is the second important property of ribbon 

operators. We will write F^ h ' 9 \t) = F^ h ' 9 \q), or, more exactly, F^ h > 9 \t) = F^ h ' 9 \q), where 
M = {x 1 , . . .,x n }. 

Linear combination of the operators F^ h ' 9 ^(t) are also called ribbon operators. They form 
an algebra T{t) = T . The multiplication rules follows 

F™(t)F n (t) = ArF k (t) Ajjy 3 = 5 hlh2 , h 5 gug 5 gM (25) 

(summation over m and n is implied). 

Any ribbon operator on a long ribbon t = t\t2 (see figure below) can be represented in terms 
of ribbon operators corresponding to its parts, t\ and t 2 



F k (ht 2 ) = n k mn F m ( tl )F n (t 2 ) n& (fcaja) = ^^h 91 (26) 

(Note that F m (ti) and F n (t 2 ) commute because the ribbons t% and t 2 do ton overlap). By 
some miracle, the tensor Q*^ is the same as in eg. (p!8[). From the mathematical point of view, 
eq. ( p6|) defines a linear mapping A(ti,t 2 ) : ^(ti,t 2 ) — > J~(ti) <%) J~(t 2 ), or just A : JF — > JF. 
Such a mapping is called a comultiplication. 

The comultiplication rules (|26|) allow to give another definition of ribbon operators which 
is nicer than eq. (f2J). Note that a ribbon consists of triangles of two types (see fig. |7|). Each 
triangle corresponds to one edge. More exactly, a triangle with two dotted sides and one dashed 
side corresponds to a combination of an edge and its endpoint, say, i and r. Similarly, a triangle 
with a solid side corresponds to a combination of an edge and one of the adjacent faces, say, j 
and /. Each triangle can be considered as a short ribbon. The corresponding ribbon operators 
are 

F^ 9 \i,r) = 6 9!l L h (i,r) F^ 9 \j,l) = T 9 ~\j,l) 

The ribbon operators on a long ribbon can be constructed from these ones. 

It has been already mentioned that the multiplication in T> and the comultiplication in T 
are defined by the same tensor Q*^. Actually, T> and T are Hopf algebras dual to each other. 
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(For general account on Hopf algebras, see |2I, [22], Q). The multiplication in T corresponds 
to a comultiplication in T> defined as follows 

A(D k ) = A™ n D m (g D n (27) 

(More explicitly, A(D {h:g) ) = J2hih 2 =h D (hi,g)® D (h 2 , g ) )• The unit element of T is \ T = £ fc F fc , 
where are given by eq. (plD ; the tensor also defines a counit of T> (i.e. the mapping 
£ : I? — > C : £(Dk) = £ k ). The unit of and the counit of are given by 

e {h ' g) = 5 g>1 (28) 

The Hopf algebra structure also includes an antipode, i.e. a mapping S : T> —>■ T> : S(Dk) = 
S^D m , or S : T -> J 7 : S(F m ) = S£\F fe . The tensor S* is given by the equation 

Sf^'f J = S a -i h h -i 5 a a -i (29) 
Here is the complete list of Hopf algebra axioms. 

Af» A*» = Ag Af" e 4 Ar = A^ = <C (30) 

nL = n* ^ mn e* n*. = e* = ^ (31) 



Af* «L = *4 fi ™ A t r A £ ^ fi L = A 1 - e* = e'e™ (32) 

S? Ap™ 1 C = # A£» = £ p e« (33) 

Most of these identities correspond to physically obvious properties of ribbon operators. 
Eq. ( |30"D is a statement of the usual multiplication axioms in the algebra namely, (F l F m )F n = 
F\F m F n ) and \F m = F m l = F m . The first equation in ([31]) (coassociativity of the comulti- 
plication in J 7 ) can be proven by expanding F k (tit2h) as Q k n F l (tit 2 ) F n (t 3 ) or Q k j F l (ti) F j (t2h) 

- the result must be the same.0 Eqs. ( |32"D mean that the multiplication and comultiplication 
are consistent with each other. To prove the first equation in (|3"2"[), expand F l (t]i,2) F m {t]t2) 
in two different ways. The second equation follows from the fact that £ q F q (tit 2 ) is the identity 
operator. 

The antipode axiom (|3~3"D does not have explicit physical meaning. Mathematically, it is 
a definition of the tensor S* : the element 7 = Sj, F k (g) Di G T ®T) is the inverse to the 
canonical element 5 = F l ® Di. The antipode have the following properties which can be 
derived from (|30|-|33|) 



S* Si A£» = Aj* S p S? n« = S" 1 S< (34) 



18 The coassociativity is necessary and sufficient for that. The sufficiency is rather obvious; the necessity 
follows from the fact that the mapping T — > is injective, i.e. the operators F k (t) with different k are 

linearly independent. 
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a) 




b) 



Figure 8: Two ribbons attached to the same site 
Finally, we can define a so-called skew antipode S* as follows 

om oi _ om oj _ K m /qir\ 
°i °n ~ D j D n — u n \°° ) 

In our case, S™ = S™, but this is not true for a generic Hopf algebra. The skew antipode have 
the following properties similar to (|33|) and (j34|) 

Sr K m Q L §k m = K Ap m n L ^ = £ P e q (36) 

~Sl Ap™ = A* S p S£ £2« = ^ S™ S< (37) 

(Note the distinction between (|33"D and fl36|), however). 

The reader may be overwhelmed by a number of formal things, so let us summarize what we 
know by now. We have defined two algebras, T> and J 7 , and their actions on the Hilbert space 
A/". In this context, we denote them by T>(a) and T{t) because the actions depend on the site 
a or on the ribbon t, respectively. Operators from T>(a) affect one particle whereas operators 
from !F(t) affect two particles. The action of T{t) on the space of n-particle states C(xi, . . . , x n ) 
depends only on the topological class of the ribbon t. This space have not been found yet, even 
for n — 2. (It will be found after we learn more about local and ribbon operators). The algebra 
T is a Hopf algebra. The comultiplication allows to make up a long ribbon from parts. There 
is a formal duality between T and T>. The comultiplication in T is dual to the multiplication 
in T>. The multiplication in T is dual to a comultiplication in T>. (The meaning of the latter is 
not clear yet). 

5.3 Further properties of local and ribbon operators 

Let us study commutation relations between ribbon operators. Consider two ribbons attached 
to the same site, as shown in fig. || a or b. Then 

F^Xu) F^ u \ qi ) = F^iqt) F( v ~ lhv > v ~ lg \ti) 

F^ 9 \t 2 ) F ( - V ' u \q 2 ) = F { - V ' u \q 2 ) F^ 9U ' lvu \t 2 ) 

In a tensor form, these equations read as follows 

F m (tO = R ik Qfj Q]£ FHqi) F l (ti) 

(38) 

F m (t 2 ) F n {q 2 ) = F f (g 2 ) F k (t 2 ) Q£ Q% R jl 
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where 

R (M(*,«) = 6h u SgA r(MM) = ^ 5 s l (39) 

Note that 

R ik fi£ R jl = R ik n% R jl = e n e m (40) 

To proved (and to see the physical meaning of) this equation, consider the configuration 
shown in fig. |^b. Clearly, F r (t 2 ti) and F s (q') commute. On the other hand, F s (q') = F s (q2qi), 
so F r (t 2 t 1 ) and F s (q 2 q 1 ) commute. It follows that R ik tt^tt^R 31 = e n e m . This identity 
can be easily written in an invariant form, namely, RR = \v®v, where R = R jl Dj ® D t 
and R = R tk Di ® D^.. It also implies that RR = lv®v because the algebra T> Cg> T> is finite 
dimensional. Thus, R = i? -1 . 

The tensor R* (or the element R € X> ® P) is called the R-matrix. It satisfies the following 
axioms 

Ajf R fcr " = R a R jn QZ R mfc A^ = fi^ R li R nj (41) 

a? = nj mr ni ns ajt R rs (42) 

where f2j mr = fi^fi^ = QY m Q l ur . Eqs. ( f4l|) follow from (|38|). Conversely, these equations 
ensure that the commutation relation are consistent. To prove the first equation in (fll|), 
commute F m (ti) F i (q 1 ) F 3 (q{) in two different ways. You will get W i3 ^ F a (q 1 ) F b (tx), with 
two different expressions for W^ b m . Then calculate W^ b m e a e b using the axioms (|30H 32^) . The 
second equation in ( f41|) can be proven in a similar way. 

To prove eq. (0), consider the configuration shown in fig. ||a. Clearly, F l (qiq 2 ) = F t (t 1 t 2 ), 
so F 3 \tit 2 ) F l (qiq 2 ) = A 3 * F k (t 1 t 2 ). On the other hand, we can first expand F 3 {tit 2 ) and 
F x (qiq 2 ) using the comultiplication rules, and then apply the commutation relations (|38|). The 
result must be the same. 

Let t be a ribbon connecting sites a and b. The local and ribbon operators commute as 
follows 

a ■ t ■ b 



F m (t) Dija) = Af gg Djja) F l (t) £>«(&) F m (t) = n£ A? F'(t) (43) 

19 This proof is not rigorous, but an interested reader can easily fix it. Anyway, you can just substitute ( |39| ) 
into (|4^) and check it directly. 



20 



These commutation relations can be also written in the form 



Dj(a)F l (t) = A? S£ n l nm F™(t) D 4 (o) 

(44) 

FtyDtf) = n l mn S£ A¥ Di {b) F m (t) 



where S* is the skew antipode (see eqs. ( p5] , |36"|) ). 

Finally, we introduce some special elements C G T> and T G T. The first one has a clear 
physical meaning: the corresponding operator C{a) = A(a) B(a) projects out states with no 
particle at the site a. The element C can be represented in the form 

C = tfDi where C (/l ' s) = iV^c^i (45) 

It has the following properties: 

CX = XC = £{X) C for any X G V S{C) = 1 

or, in tensor notations, 

ft* c* = fijjc* = &,-c* e fc c fc = 1 (46) 

The element T G is dual to C; it is defined as follows 

T = Ti F l where r (M = iV" 1 ^ (47) 

Its properties are as follows 

A^'r, = A?Ti = e J 'r fc e fc r fc = 1 (48) 

Note that Tfc C fc = A^~ 2 . Using these properties, we can we can derive an important consequence 
from the commutation relations (|43|) 

T s {l s mp S p q F m (t)C(a)Fi(t) = T s n s pm S p q F«(t)C(b)F m (t) = N~ 2 (49) 



5.4 The space C(a, b) 

Now we are in a position to find the space C(a, b) and to prove the assertions from Sec. |5.1| . 
The first assertion was that any excited spot can be transformed into one particle. It is simple 
if we can transform two particles into one by ribbon operators. Let us choose an arbitrary site 
b the excited spot to be compressed to. Let some constraint, A(s) — 1 = or B(p) — 1 = 0, 
be violated. Choose any site a containing the vertex s or the face p. Connect a and b by 
a ribbon. By the assumption, we can clean up the site a while changing the state of b, but 
without violating any more constraint. We can repeat this procedure again and again to clean 
up the whole spot. 

So, we only need to show that two particles can be transformed into one. What does it mean 
exactly? Physically, any transformation must be unitary, but it can involve also some external 
system. (Otherwise, it is impossible to "decrease entropy", i.e. to convert many states into 
fewer). On the other hand, it is clear that unitarity is not relevant to this problem. However, 
we should exclude degenerate transformations, such as multiplication by the zero operator. So, 
it is better to reformulate the assertion as follows: any two-particle state (plus some other 
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excitations far away) can be obtained from one-particle states (plus the same excitations far 
away). Let \ip) G £(a, b, . . .) be such a two-particle state. We are going to use the formula fl49[). 
Let 

G q = N 2 r s n s mp S p q F m (t) \ v q ) = C{a)F\t)\^) (50) 

Then = G q \i] q ). The states \rj q ) belong to £(b, . . .), i.e. do not contain excitation at a. 
This is exactly what we need. 

The other two assertions were about the action of local operators on the space £(a,b), so 
we need to find this space first. We can consider this space as a representation of the algebra 
£ = £{t) generated by the operators Dj = Dj(a), F l = F l (t) and Dj = Dj(b). As a linear 
space, £ = T> <g> JF <g> V. (Thus, £ has dimensionality N 6 ). Multiplication in £ is defined by the 
commutation relations ([|3J). We will call £ = £(t) the algebra of extended ribbon operators. It 
is just an algebra, not a Hopf algebra. More exactly, it is a finite-dimensional C*-algebra. The 
involution (=Hermitian conjugation) is given by the formulas (cf. (|19|)) 

(/W f = JVW-*) (^ (M ) f = F [h ~^ 9) iD{ h J = ( 51 ) 

[Remark. Apparently, the algebra £ will play the central role in a general theory of topolog- 
ical quantum order. Indeed, we were lucky to define ribbon operators separately from local 
operators. In the general case, ribbon operators should be mixed with local operators.] 

So, we are looking for a particular representation £ of the algebra £. This representation 
must contain a special vector |£) (the ground state) such that 

D k \0 = e k \0 D' k \0 = £ k \0 (52) 

We start with constructing a representation £ spanned by the vectors \ip k ) = F k (It will be 
proven after that C = £). We assume that the vectors \%l) k ) are linearly independent. This need 
not be the case in the representation £ but we can postulate \ip k ) being linearly independent 
in L. Thus, £ contains a factor-representation of £. 
The representation £ is given by the formulas 

DjW) = S^JUVO Fi\^ k ) = A^|V m > D'^ k ) = n^liT) (53) 

It is easy to show that this representation is irreducible. Hence, £ contains £, i.e. the vectors 
\%l) k ) are linearly independent in £. The scalar products between the vectors \ip k ) can be found 
from (|53|) and QSiP , 

^(^^(M) = N- l 8 v>h 6 u>g (54) 

To prove that £ = £, we use the equation (|49| ) again. For an arbitrary two-particle state 
\ip) G £, define the vectors \rj q ) and operators G q as in eq. (|50|). Then \tp) = G q \r] q ). One 
could say that \r] q ) G £{b) but, actually, the space £{b) is spanned by the sole vector It 
follows that \tp) G £ — the assertion has been proven. Thus, the action of local and ribbon 
operators on the space £ = £(a, b) is given by eq. (|53|) . 

It is easy to see that the action of T>{a) on £(a, b) is exact (though it is reducible). Besides 
that, T>(a) is the commutant of T>(b) in L(£(a, b)) and vise versa. (That is, T>(a) consists of all 
operators X G L(£(a, b)) which commute with every Y G T>{b) ). Hence, T>(a) includes all local 
operators acting on the space £(a,b). Indeed, a local operator, which involves only spins near 
the site a, must commute with any operator acting on distant spins. Of course, there are many 
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such operators, but their action on the two-particle space £(a, b) coincides with the action of 
the operators from T>(a). This is also true for a multi-particle excitation space C(xi, . . . ,x n ). 

The space C{x\, . . . , x n ) can be described as follows. Let us connect the sites Xx, . . . , x n by 
ri — 1 ribbons ti, . . . , t n _\ in an arbitrary way so that the ribbons form a tree. Then the vectors 
\ip kl '---' k n~i} = F kl (ti) . . . F fe "- 1 (t n _ 1 ) |£) form a basis of C(x 1 , . . . , x n ). Choosing different 
ribbons means choosing a different basis. In the next section we will give another description 
of multi-particle excitation spaces. 

6 Topological operators, braiding and fusion 

Let us consider again the n-particle excitation space C = C(x\, . . . , x n ). The algebra L(£) 
includes the local operator algebras T>(xi), . . . ,T>(xi). An operator Y G L(£) which commute 
with every X G D(xj) (j = l,...,n) is called a topological operator. Physically, topological 
operators correspond to nonlocal degrees of freedom. For n = 2, the algebra of topological 
operators coincides with the center of T>(x\) or T>(x2). (The two centers coincide). Hence, the 
only nonlocal degree of freedom is the type of either particle. (The two particles correspond to 
dual representations of T>\ in other words, these are a particle and an anti-particle). So, there 
is no hidden (i.e. quantum nonlocal) degree of freedom in this case. Such hidden degrees of 
freedom appear for n > 3. 

To describe the space C and operators acting on it, let us choose an arbitrary site Xq 
(distinct from x±, . . . , x n ) and connect it with xi, . . . , x n by non-intersecting ribbons ti, . . . , t n , 
see fig. |lOla. As stated above, the space C(xq, x\, . . . , x n ) is spanned by the vectors 

|^ fcl '---' fc ") = F fcl (ti)...F fc "(tn)|0 (55) 

The space in question, £ = C(xi, . . . ,x n ) is contained in C(xo,x\, . . . ,x n ). It consists of all 
vectors \ip) G C(xo, Xi, . . . , x n ) which are invariant under the action of V(xq) on the latter 
space. 

The advantage of this description is that we can easily find all operators on the space 
£(xo, x\, . . . , x n ) which commute with V(x\) ® . . . ® V(x n ). These are simply operators which 
act on the ends of the ribbons ti, . . . , t n attached to the site Xq. More exactly, an operator D^p 
(r = 1, . . . ,n) acts on the r-th ribbon as Dj = Dj(xo) (see eq.(|5B|)), but does not affect the 
other ribbons, 

Thus we arrive to an interesting physical conclusion. Let us consider only one particle attached 
to an end of a semi- infinite ribbon (an analog of Dirac's string). Then the topological operators 
act on the far end of the ribbon. 

Example. Let us see how the topological operators act on magnetic vortices. As shown in 
Sec. pTT| , a vortex type is characterized by a conjugacy class C of the group G. Individual 
topological states of the particle are characterized by particular elements v G C. In terms of 
the notation (|53D, such a state can be represented as follows 

\ U)V ) = |ci i/2 \^ (u,x) ) 

x: x —1 ux=v 
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a) b) c) 

Figure 10: Braiding and fusion in terms of ribbon transformations 

where u G C characterize the local state of the particle. One can easily check that D', h g , \u, v) = 
8h,gvg-A u i This is consistent with eq. (p2j). Note that the local degree of freedom, it, is not 
affected. 

How can we physically apply topological operators to particles? We can just move the 
particles around each other; this process is called braiding. Let us see what happens if we 
interchange two particles, x s and x s+ \, counterclockwise, as shown in fig. [To|b. The state 
1^,. . . , fc, i, . . \) ecomes a new state 

To represent this state in the old basis, we should represent the operator F k (t' s ) F l (t' s+1 ) in 
terms of F m (t s ) and F n (t s+1 ). Obviously, F k (t' s ) = F k (t s+1 ); also F l (t' s+1 ) = F l (t s ) as long as 
there is no particle at x s+ i, i.e. the operator F fc (t s+1 ) is not applied yet. Hence 

F k (Q F\t' s+1 ) = F k (t s+1 ) F l (t s ) 

Now we can apply the second commutation relation from (p8[). (Actually, we should reverse 
it). It follows that 

F k (t s+1 ) F\t s ) = R ji n l mi n k j F m (t s ) F n (t s+1 ) 

(see eq. (|56|)). Consequently, the counterclockwise interchange operator has the form 

K„ = R ji (D'i <g> £>;.) a = a R ij (D^ ® Dj) (57) 

where G is the permutation operator, and D' i: Dj are understood as topological operators. 
(Note that the operator o permutes both topological and local degrees of freedom). 

Example. Consider two magnetic vortices characterized by topological parameters v i, v 2 G G. 
The operator 1Z -r\ acts on the state \vi,v?) as follows 

K^\vi,v 2 ) = \viv 2 v^, vi) (58) 

(The local parameters, u\ and u%, are suppressed in this formula). 
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Finally, let us see what happens if two particles, x s and x 8 +i, fuse into one. The resulting 
particle can be characterized by the action of topological operators on it. From this point of 
view, we can just glue parts of the corresponding ribbons (see fig. pX|c) instead of fusing the 
particles themselves. Then 

F k (t s ) F l (t s+1 ) = n l nj A* F™(t' s ) F n (t' s+1 ) F*(q) 

where .D^. acts on the end of the ribbon q. Thus, fusion is described by the comultiplication in 
the algebra V, see equation ([27]). (To avoid confusion, one should replace D± with in that 
equation). The topological operator A(D' k ) acts on a particle pair as the topological operator 
D' k on the particle resulting from fusion. 

Example. Consider a pair of opposite magnetic vortices \v, t> -1 ). The operators A(D k ) act 
on this state as follows 

A(D{ htg) ) \v, v- 1 ) = <f M \gvg-\ gv^g' 1 ) (59) 

It terms of the representation classification (see Sec. |5.1|) , this action corresponds to the pair 
(C, x), where C = {1}, and x is the adjoint representation of G. Thus, when opposite magnetic 
vortices fuse, the resulting particle has no magnetic charge but may have some electric charge. 

7 Universal computation by anyons 

(This section should be considered as an abstract of results to be presented elsewhere). 

Universal quantum computation is possible in the model based on the permutation group 
S5. (Unsolvability of the group seems to be important). Vortex pairs \v, v^ 1 ), where v is a 
transposition, are used as qubits. It is possible to perform the following operations. 

1. To produce pairs with zero charge. If a pair is created from the ground state, it has no 
charge automatically. 

2. To measure the electric charge of a vortex pair destructively. For this, we should simply 
fuse the the pair into one particle. 

3. To perform the following unitary transformation on two pairs 

\u, m -1 ) <8> \v, v~ x ) 1— > \vuv _1 , vu^v^ 1 ) (g) \v, v~ x ) (60) 
For this, we pull the first pair (as a whole) between the particles of the second pair. 

4. To measure the value of v and produce an unlimited number of pure states \v, v~ l ) for 
any given transposition v (say, (1, 2) or (2, 3) ). [At first sight, it is impossible because we 
can only measure the conjugacy class of a v. However, we can agree on a given state to 
correspond to v — (1, 2). Then we use it as a reference to produce an unlimited number 
of copies.] 

The operations || and |4] are sufficient to perform universal classical computation. It is 
relatively simple to run quantum algorithms based on measurements [^] . Simulating a universal 



gate set is more subtle and requires composite qubits. That is, a usual qubit (with two distinct 
states) is represented by several vortex pairs. 
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Concluding remarks 



It has been shown that anyons can arise from a Hamiltonian with local interactions but without 
any symmetry. These anyons can be used to perform universal quantum computation. There are 
still many things to do and questions to answer. First of all, it is desirable to find other models 
with anyons which allow universal quantum computation. (The group S5 is quite unrealistic for 
physical implementation). Such models must be based on a more general algebraic structure 
rather than the quantum double of a group algebra. A general theory of anyons and topological 
quantum order is lacking. [In a sense, a general theory of anyons already exists JTU] ; it is based 
on quasi-triangular quasi-Hopf algebras. However, this theory either merely postulates the 
properties of anyons or connects them to certain field theories. This is quite unlike the theory of 
local and ribbon operators which describes both the properties of excitations and the underlying 
spin entanglement.] It is also desirable to formulate and prove some theorem about existence 
and the number of local degrees of freedom. (It seems that the local degrees of freedom are a 
sign that anyons arise from a system with no symmetry in the Hamiltonian). Finally, general 
understanding of dynamically created, or "materialized" symmetry is lacking. There one may 
find some insights for high energy physics. If we adopt a conjecture that the fundamental 
Hamiltonian or Lagrangian is not symmetric, we can probably infer some consequences about 
the particle spectrum. 
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